PhysicsO.Ol: Object-Oriented Programming for Exact Diagonalization 
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A new system of library code is proposed and initiated. It is emphasized that the same terminolo- 
gies as we find in our textbooks should be used for class names in the library code. The language 
C# invented by Microsoft is adopted in this project. Several rules of thumb are suggested in order 
to obtain easy-readable coherent codes. As a first step, we present the library code for exact diago- 
nalization in physics. When we build codes, we clearly distinguish between model independent and 
dependent parts, and we use familiar terminologies like Hamiltonian, HilbertSpace, GroundState, etc 
as class names. As an explicit example, we calculate ground state energy of a quantum dot, showing 
the triplet-singlet transition. 
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I. INTRODUCTION 

In Feb. 2002, Microsoft launched a new program lan- 
guage named C#, which is strongly object-oriented^]]. 
It is argued that C# language has three characteristic 
features, 1. C# is as elegance as Java, 2. C# is as 
powerful as C++, 3. C# is as productive as Visual Ba- 
sic. Although it is still questionable whether or not C# 
is useful for scientific computing in academic area, we 
here adopt C# because of its strong character of object- 
oriented language. 

Reusability is always one of main concerns in the area 
of software. Scientific computing also should require 
reusability seriously. We notice many efforts to improve 
reusability. The web site managed by Troyer 2J is re- 
markable. The library code for density matrix renormal- 
ization group method invented by White is also well 
known. These library codes are written in CH — h There 
are many fortran codes in netlib||. It is clear that there 
is a general trend of transition from non object-oriented 
to object-oriented language. 

The main advantage of C# over C++ is efficiency 
in making Windows applications. Since computer users 
are addicted by Windows operating system, user friendly 
programs would be Windows applications. As far as end 
users need more user friendly programs, C# can be use- 
ful. One good point is that so called Mono project is 
carried on to develop the C# compiler for Linux. After 
successful finish of Mono project, we guess that C# will 
become more popular. 

In this article, we propose a new way of making li- 
brary codes with C# as a preparation for the battle be- 
tween C# and C++, which will take place soon. At 
present, CH — h is the most powerful. However, after some- 
one makes ease access to parallel computing in loosely 
connected Windows operating systems using XML web 
services, we expect that C# can compete with C++ even 
in scientific computing. 
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II. SUPERCODE 

In order to improve the power of scientific computing, 
we have to consider the three factors: CPU time, REM 
memory, and Coding time. For some complicated prob- 
lems, the main factor is nothing but coding time. In this 
situation, someone want to use Standardized Reusable 
Components (SRC). When we say Hamiltonian, everyone 
in physics community can understand the meaning of it. 
Hence, there must be a class called Hamiltonian. At this 
moment, we should notice the usefulness, when we use 
the same terminology as we find in our textbooks. All li- 
brary codes must be organized like our textbooks and our 
library. The system of library codes should not be in sim- 
ple alphabet order, but would be in subject-based order. 
This subject-based integration of library codes will result 
in high correlations between codes, in other words, it will 
need more coherence between codes. It is not an easy task 
to make classification in each category of our knowledge. 
We may need a big discussion on this matter. At any 
rate, this library codes can be called Supercode. As we 
see in usual real library, we can divide Supercode as 
Biology, Chemistry, ElectricalEngineering, Math- 
ematics, MechanicalEngineering, Physics, etc. 

In scientific computing, programmers always want to 
know results as quickly as possible. Because of this 
hurry habit, program codes become less structural, and 
reusability is lost. In order to enhance reusability, it is 
strongly proposed that scientists should make program 
codes in distinguishing model-independent and model- 
dependent parts. The model-independent parts are prop- 
erly organized and will be upgraded in the future. Nu- 
meric workers should be familiar with the whole structure 
of model-independent parts, and should handle model- 
independent and model-dependent parts simultaneously. 
In order to emphasize the way of coding, we denote 
HANA, which is abbreviation for "Handle All aNd All" , 
where the first "All" means model-independent parts and 
the second "All" means model-dependent parts. We de- 
fine HANA project as efforts to make Supercode. We 
are summing up some rules of thumb for Supercode: 

1. Divide and Conquer. Do not exceed the total num- 
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ber of code over 100 lines for each class. 

2. Handle all and all. Consider always a general code 
for a general case as much as possible. Collect the codes 
of model-independent parts for Supercode. 

3. Give proper full names. Be careful in choos- 
ing namespace names and class names. For instance, 
all codes in Numerical Recipes should become private 
methods in Supercode. An easy way to choose names 
is to use an index in the end pages of textbooks. Allow 
no space in a full name, instead use upper case charac- 
ters, for instance, ProteinFolding, SpecialFunction, 
QuantumMechanics, etc. This rule was adopted in the 
textbook written by Deitel, "C# How to program" . 

4. Use SRC. Use confirmed codes as much as possible. 
In order to avoid any mistakes, do not use special or 
home-made codes. 

5. Reduce number of arguments in constructors and 
public methods. Future library users will definitely want 
a simple structure of arguments. 

6. Consider polymorphism in the pattern of "word 
after the same word". Make an abstract class for the 
common name. For instance, when we note Transla- 
tion Symmetry, Rotation Symmetry, Reflection Symme- 
try, Spin Inversion Symmetry, and so on, we make ab- 
stract class Symmetry and subdirectory Symmetry, and 
then make class Translation, Rotation, Reflection, Spin- 
Inversion in Symmetry directory, using the key word 
override. 

Because of the author's speciality, only a small part 
of Supercode is discussed in the following. The part is 
named as PhysicsO.Ol. The version number 0.01 repre- 
sents the first step toward the final goal of Supercode. 
We hope that upgrade version of Physics will be used 
for all computing processes in physics. 



III. PHYSICS0.01 



We here omit to show the present codes of 
PhysicsO.Ol. All detailed codes are freely open in the 
author's home pageQ- This prototype object-oriented 
programming for physics is invented mainly only for 
exact diagonalization. We expect that Monte Carlo 
and Maxwell equation will be included. Hence upgrade 
Physics0.02 will appear soon. 



IV. EXAMPLE: GROUND STATE ENERGY IN 
2D PARABOLIC QUANTUM DOT 

In order to present the general procedure of HANA 
project, we consider the problem of evaluating ground 
state energy in a two-dimensional parabolic quantum dot. 
In this problem, the second quantized Hamiltonian[7J is 
written as 

H = ^2 {hjj{2n + \l\ + 1) - -hjJ C l ~ g^haj c o-)cl li:r c n i a 

7i, I. a 
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where one-particle creation operators c\ la have three 
quantum numbers^. The principal quantum number 
n runs as 0, 1, 2, 3, • • • , while the angular momentum 
quantum number I is given by 0, ±1, ±2, ±3, ■ ■ ■ , and 
the spin index a = ±1. Note that the confining potential 

is related to ujq in the usual way, and u — \Jlo$ + \^>c, 
where the so called cyclotron frequency is proportional 
to the strength of magnetic field B as ujq — eB/m*c. It 
is straightforward to calculate the Coulomb matrix ele- 
ments: 



Following the above rules, we try to construct 
PhysicsO.Ol. All files related to physics exist in a single 
directory named PhysicsO.Ol, which contains many sub- 
directories and files. Here the subdirectories are names 
of namespaces, and files are names of classes. The sub- 
directories, which can be called volumes, are classified 
as usual as textbooks for physics. We make the sub- 
directories: ClassicalMechanics, QuantumMechan- 
ics, Electromagnetism, StatisticalMechanics, Gen- 
eralPhysics, etc. One of these volumes, Quantum- 
Mechanics, has many class files: ConservedQuantity.es, 
GroundState.es, etc, and subdirectories: Hamiltonian, 
HilbertSpace, Symmetry, Operator. The directory 
Hamiltonian contains many classes. Following the sixth 
rule of polymorphism, the class files in Hamiltonian 
directory are given by FermionOneBody.es, FermionT- 
woBody.cs, Heisenberg.es, Kondo.cs, and abstract class 
Hamiltonian.es, since we note Fermion One Body Hamil- 
tonian, Fermion Two Body Hamiltonian, Heisenberg 
Hamiltonian, and Kondo Hamiltonian. Still many other 
class codes should be made, accumulated and organized 
in QuantumMechanics of PhysicsO.Ol. 



1 1 6 ui*lu 

2 ^ / niri2n3«4 (^) ^) m ) — 2~K \ 2H ^"i™ 2 ™ 3 ™ 4 ^' m ^ 



n \ m 



58.47 m* lmeV 
nio \l — \l — C nin2n3ni (k,l;m), (2) 



where k is a dielectric constant, and C ni „ 2ri 3 ri4 (fc, I; m) 
are dimensionless numbers. 

With the usual bulk values, m* — 0.07m and k = 13, 
the overall factor of the Coulomb interaction becomes 
58.47/13 * V0.07 = 1.19. However, since the used quan- 
tum dot looks two-dimensional, it seems that the values 
of dielectric constant and effective mass are not known. 
Thus, it is needed to take the overall factor as an input 
parameter. The value of loq is also an input parameter. 
In fact, the side gate is used in order to control the con- 
fining potential 9]. In consequence, this calculation has 
four input parameters: number of electrons, the overall 
factor, cuq, and cue- We take 20 lowest energy states as 
2n + \l\ + 1 < 4 in this exact diagonalization. Only the 
main parts of the code arc presented as follows. 
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us ing QM=Phy s i c s . QuantumMe chani c s ; 
using H=Physics . QuantumMechanics . Hamiltonian; 
using HS=Physics . QuantumMechanics . HilbertSpace ; 
using OP=Physics . QuantumMechanics . Operator; 



int Ne = Convert . ToInt32 (textBox2 . Text) ; \\ input number of electrons 
double factor = Convert .ToDouble(textBox3 .Text) ; \\ input the overall factor 
double omegaO = Convert .ToDouble (textBox4. Text) ; \\ input omegaO 
double omegac = Convert .ToDouble (textBox5 .Text) ; \\ input omegac 

int half = (numberOf Electrons + (numberOf Electrons °/„ 2))/2; \\ At half, S_{z}=0 

QM. Space il = new QM.Space(7,2) ; \\ 14 sites 
QM. Space i2 = new QM.Space(3,2) ; \\ 6 sites 

QM. Space ii = new QM.Space(il,i2) ; \\add Space il and i2, so 20 states 

QM . QuantumNumber qn = new QuantumNumberForModel(ii) ; Wassign n, 1 and spin 
QM. Coupling u = new CouplingForModelO (qn, omegaO, omegac) ; \\one body spectrum 
QM. Coupling w = new CouplingForModell (qn, factor, omegaO, omegac) ; WCoulomb matrix 

OP . FermionCreation cdagger = new OP . FermionCreation(ii) ; 
OP . FermionAnnihilation c = new OP .FermionAnnihilation(ii) ; 

H . Hamiltonian hO = new H.FermionOneBody(u, cdagger, c) ; 

H . Hamiltonian hi = new H.FermionTwoBody(w, cdagger, cdagger, c, c) ; 

H.TotalHamiltonian h = new H . TotalHamiltonian(hO ,hl) ; 
h. ZeroPointEnergy = 0.0; 

OP. Number nop = new OP. Number (numberOf Electrons) ; 

for(int ta = 0; ta < 3; ta++){ 

for (int ts = half; ts < half +2; ts++){ 

QM.ConservedQuantity totalAngularMomentum = new QM. ConservedQuantity (qn, 1 ,ta) ; 
QM.ConservedQuantity totalSpin = new QM. ConservedQuantity (qn, 2, ts) ; 

HS . HilbertSpace hill = new HS . Fermion( cdagger , nop); 

HS . HilbertSpace hil2 = new HS . Reduction (hill .totalAngularMomentum) ; 

HS . HilbertSpace hil = new HS. Reduction(hil2, totalSpin) ; 

QM.GroundState psi = new QM.GroundState(h,hil) ; 
textBoxl . Text = Convert .ToStr ing (psi .Energy) ; }} 



Note that the classes of model-independent parts 
are in specific directories. The model-dependent parts 
consist of QuantumDot.es, QuantumNumberForModel.es, 
CouplingForModelO.es and CouplingForModell.es. The 

user friendly simulation program runs in Windows op- 
erating system. Giving the four input parameters, we 
simply click a button to find the ground state energy. 
The final result of ground state energy is printed on the 
screen. 

We present the result of calculation related to the 
triplet-singlet transition|lOj. when the strength of confin- 



ing potential is increased in the dot containing six elec- 
trons. We let the overall dimensionless factor be 2, and 
uic = 0. We find the transition near Hloo = 1.14 meV as 
shown in Table [I] We expect that someone observe this 
transition using the side gate voltage in the future. 

V. CONCLUSION 

We proposed a library code system, emphasizing that 
class names should be the same as we find in textbooks. 
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TABLE I: The ground state energy in each sector. We notice 
that there is the triplet-singlet transition between 1.1 and 1.2. 



Sector 




Ground State Enerj 


?y 


Sz Lz 


huo — l.i 


huio = 1.15 


huuo — 1.2 





59.7012 


61.4422 


63.1693 


1 


59.7247 


61.4975 


63.2490 


1 


59.7012 


61.4759 


63.2299 


1 1 


59.7247 


61.4975 


63.2490 



Like an well known idiom in algorithm, "Divide and con- 
quer" , we introduced "Handle all and all" for this library 
system. The first "all" means library codes and the sec- 
ond "all" means our specific project codes. Hence, there 
are only two folders: one is library and the other is our 
specific model dependent project. As the library grows, 
we expect that each class becomes strongly correlated 
with others. To change codes is not an easy task because 
of coherence between codes. However, experts can freely 
change the codes, and thus invent the upgrade version 
with fixing the first number such as Physicsl.03 from 



the mother version Physicsl.OO. For instance, if some- 
one feel that the class Symmetry can be improved, then 
change the code. If this change is worthwhile for ev- 
eryone, then it should be reflected in the next version 
of Physics with clear notification of contributors. Along 
this line, the prototype library code, PhysicsO.Ol, will 
be upgraded soon. 

We know that nowadays many theoretical physicists 
spend more time in complicated computing, which is 
however conceptually trivial. We agree that physicists 
need to do philosophical and conceptual thinking more. 
The author likes to imagine that a higher version of 
Physics will reduce the time of computing for all physi- 
cists. 
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